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Abstract 

Context. Magnetic element tracking is often used to study the transport and diffusion of the magnetic field on the solar 
photosphere. From the analysis of the displacement spectrum of these tracers, it has been recently agreed that a regime of 
super-diffusivity dominates the solar surface. Quite habitually this result is discussed in the framework of fully developed 
turbulence. 

Aims. But the debate whether the super-diffusivity is generated by a turbulent dispersion process, by the advection due to the 
convective pattern, or by even another process, is still open, as is the question about the amount of diffusivity at the scales 
relevant to the local dynamo process. 

Methods. To understand how such peculiar diffusion in the solar atmosphere takes places, we compared the results from two 
different data-sets (ground-based and space-borne) and developed a simulation of passive tracers advection by the deformation 
of a Voronoi network. 

Results. The displacement spectra of the magnetic elements obtained by the data-sets are consistent in retrieving a super- 
diffusive regime for the solar photosphere, but the simulation also shows a super-diffusive displacement spectrum: its competitive 
advection process can reproduce the signature of super-diffusion. 

Conclusions. Therefore, it is not necessary to hypothesize a totally developed turbulence regime to explain the motion of the 
magnetic elements on the solar surface. 
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1. Introduction 

The outermost shell of the Sun is dominated by turbulent 
convection, which is the main driver of heat and entropy 
exchange between the different layers. The present 
analytical models of convective turbulence (namely, 
the Mixi ng Length Theory and the Full Spectrum 
Theory - Bohm -Vitense 1958 : Canuto fc Mazzitellilll99ll 


(Romano et al. 2012: Kevs et al. 

2014) and they (or their 

proxies, see: ISteiner et al.l 2001: 

IJitenbroek & Tritschler 


20061) are often used as tracers of the plasma flows (e.g., 
2008t Berrilli et al. 2013 . 2014). Moreover, 


Bonet et al 


Canuto et al.l 1996h arenot able to describe appropriately 
the convective shell from the small subgranular scales 
(< 100 km) to the global scale . On the other ha nd, the 
fully numerical aproaches (e.g., Beeck et al.ll20ld) . while 
succesfull in reproducing limited regions of the Sun, 
are also unable to provide a complete view, due to the 
enormous range of temporal and spatial scales needed 
to simulate the whole solar convective envelope and the 
limited computing power available. 

On the observational side, a viable way to study the 
convective flows on the solar surface are the small 


the study of the dynamics of the magnetic elements on the 
solar photosphere is of particular interest in astrophysics, 
since it could also provide constraints on: 

a) the characteristic spatiotemporal scales of the emer- 
gence and diffusion proce s ses of the magn etic field (e.g., 
lOrozco Suarez et"^ 12012 : Stangaliiiil 2014 ) 

b) the rates of interaction between the magnetic elements 
themselves, which can cause magnetic reconnections and 
subsequent micro-flaring events: a plausible mechanism 
for the heating of the u pper solar atmosphe re (as, for 
example, in the model of IViticchie et al. 20061) . 


Send offprint requests to: delmoro@roma2.infn.it 


In order to simplify the description of the motion of 
the magnetic elements, we neglect the force they exert 
on the surrounding plasma ( Petrova^ 1994 . 2001 ). thus 
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considering them as passive tracers of the underlying 
flows. This assu m ption ha s been verified in the quiet Sun 
dViticchie et al. _ 20101: Orozco Suarez fc Bellot Rubiol 


Bellot Rubio k, Orozco Snared 

2 OI 3 II . where high 


ioii 

values of /3 

are observed. 

In a Lagrangian approach, the trajectories of single 
magnetic elements are employed for a statistical deter¬ 
mination of the dynamical properties of the ensemble. 
In particular, the magnetic element motion is described 
via the mean square displacement and modeled with 
a power law: {5"^) = k ■ V, where (6^) is the mean 
square displacement of the magnetic elements from their 
initial (t = 0 ) position, k plays the role of a diffusion 
coefficient, t is the time span from the magnetic element 
hrst detection, and 7 is the so-called spectral index of the 
power law. For diffusive motions, if 7 =1, the diffusion 
is termed normal or Fickia n, (as, e.g., in the case of the 


20121 : 


Giannattasio et al. 


Random Walk scenario - lEinsteinl 11905!) . When 7 ^1 


the diffusion is termed anomalous. Sub-diffusion if 7 <1 
(the diffusion process is “slower” than normal diffusion); 
super-diffusion if 7 >1 (the process is “more rapid” than 
normal diffusion). It is worth to note that only in the 
7 = 1 case the constant k coincides with the diffusion 
coefficient. In the anomalous diffusion cases, the diffusion 
coefficient, defined as k{Sr,t) = ( S^)/t, depends on the 


spatial and temp (3 ral scale (e.g., ISchriiyer et al. 19961 


Chae et al 

2m2li . 


2008 : Abramenko et al. I I 2 OII : Lenreti et al 


The diffusion of magnetic elements has been studied 
by several authors through segment ation and tracking 


2012 . 


algorithms of photosperic data (see iGoode et al 
and references therein). 

The magnetic elements are detected using either 


the associated polarization (e.g., IViticchie et al, 
2 OO 9 II or their excess brigh tness in the G-band (e.g., 
Sanchez Almeida et al. 2010l) : then, most authors utilize 
a Lagrangian approach to determine t he spectral index 


7 (see, e.g., Ruzmaikin et al. lOOfil : IJtz et all 2010l: 


Kitiashvili et al. 2012ll . See for example Jafarzadeh et al.l 
(I 2 OI 4 for the most recent review on the results obtained 
in previous works. 

Apar t from few older reports (lSchiTiver_^Mar^ 


199C; Lawrence fc Schriiver 199,4 Gadavid et al. 19981 


199911 , all the works 


Sion 


on 


on magnetic 
g a super-dii 


elements diffu- 
hsive behaviour 


(e.g. 

Wang & Sheelev 

1993: Schriiver et al.l 

1994 

BerE 

'■er et al. 1998: Hagenaar et al. 1999: Ghae et al. 

2008 

c Abramenko et al. 

201ll Ghitta et al. 

2012 : 

Giannattasio et al. 201,4 

2014a: .lafarzadeh et al. 

20141 


at all accessible spatio-temporal scales, ranging from 
the granu lation (diam e ter cs IMm and lifetime cs few 
minu tes; lMulleilll999l : Hirzberger et al. 19991 : Del Morol 
2004 to the supe r-granulation (diameter 30Mm and 


lifetime a day; Raiu et al. 19991 : Srikanth et al. 20001: 


Del Moro et al. 2 OO 4 I 1 . 


In the literature (e.g., Abramenko et al.l 201lh . ((5^) 


is 


often called “average squared displacement” or “squared 


displacement spectrum” or simply “displacement spec¬ 
trum”. In the following, for the sake of clarity and to 
be consistent with previous literature we will use the 
term displacement spectrum when referring to the (5^) 
computed from the magnetic element tracking. 

In this work, we used an advection simulation to repro¬ 
duce the scaling laws observed in the magnetic element 
motion on the solar photsphere. We assume that the 
plasma velocity structures (i.e.: the granules) evolve under 
a simple mutual repulsion model, and that the dynamics 
on the granular scale is completely determined by the 
deformation of such granular structures. We consider 
the connected down-flows on the solar photosphere as a 
sort of lattice-like network, and suppose the magnetic 
elements to be bonded to the photospheric velocity sinks 
located at the lattice junctions (the intergranular lanes 
vertexes). Under such hypotheses, the simulation returns 
a “super-diffusive” regime with a scaling similar to that 
observed on the solar surface. 

Therefore, a signal of super-diffusion can be imitated as¬ 
suming that the magnetic elements are passively advected 
by the plasma down-flow sinks and that the evolution 
of the granular pattern is a competitive expansion process. 


2. Real data 


In Figs. [T] and [ 2 j we report the displacement spectrum 
(Sf) for two different data-sets imaging the quiet Sun. 
The plot in Fig. [1] has been realized computing the 
displacement of G-band bright points co-spatial with 
magnetic field, as imaged by the ground based instrument 
IB IS (Interferome t ric B l dimensional Sp e ctropolarimeter 
Gavallini et al. 200 it Gavallinil 2004 Viticchie et al. 


2009ll . located in the Dunn Solar Telescope of the National 


Solar Obs ervatory (NM-USA). Th is analysis is similar 


to that in lAbramenko et al.l (l2011ri , the main difference 


being the check of the magnetic origin of the G-band 
bright points allowed by the co-temporal magnetograms. 
The spatial resolution of this data-set, restored 
with MFBD (Multi-Frame Blind Dec onvolution 


2004 . 


IS 


Van Noort. Rounne van der Voort fc Lofdahl 
~ 0 ” . 1 , and the time between consecutive frames is dt = 5 
s. The data have been kh — filtered with a 7 km/s 
cut-off velocity. Since data are affected by seeing, we tried 
to compensate this effect by applying a jitter removal 
algorithm. We subtracted from each bright point dis¬ 
placement, a tip-tilt contribution. Such contribution was 
estimated via the correlation of a 3” .3 x 3” .3 sub-frame 
centered on the bright point at two consecutive frames. 
The plot in Fig. [2] has been realized computing the 
displacement of magnetic elements from the long- 
duration (~ 25 hr ) data-set acqu ired by HINODE-NFI 
( Kosugj et al.1 200^ iTsuneta et al.l 2008h an d used in the 
analysis bv iGiannattasio et al. ( 201,4 2014a[) . The spatial 
resolution of this data-set is ~ 0”.3 and the time between 
consecutive frames is dt = 90 s, also these data have been 
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Figure 1. The displacement spectrum for magnetic G- 
band bright points for a high-resolution data-set acquired 
at the NSO-DST observatory. The lines represent power 
law fits for the ranges [t < 40 s] s and [40 s < t < 140 s] . 


Figure 2. The displacement spectrum for magnetic ele¬ 
ments from a long-duration spectropolarimetric data-set 
acquired by HINODE-NFI. The lines represent power law 
fits for the ranges [90 s < t < 1600 s] and [1600 s < t < 
10000 s]. 


kh — uj filtered (f or further details on the calibration see 
Gosic et al ][ 20 il). 


The tracking alg orithm used for th e two data-se ts is 
the same us ed in Giannat tasio et al. (I 2 OI 3 I l2014al) a nd 
presented in lPel Moro (2004) and Berrilli et al.l ( 20051) 
These two data-sets combined allow us to explore more 
consistently the time range from 5 s to ~ 10000 s. They 
consistently identify a super-diffusive behaviour (7 ~ 1.4) 
in the overlapping time range 90 s < t < 500 s. We stress 
that, to our knowledge, this is the first time that two 
data-sets acquired with rather different instruments and 
at different sites, agree on retrieving 7 values consistent 
within the errors. 


3. Data Interpretation 

On the light of these two plots, whose results are summa¬ 
rized in Table HI we put forward a possible interpretation. 
For times shorter than 40 s (and hence for space scales 
smaller than ^ 50 km), the displacement behaviour looks 
slightly super-diffusive: 7 ~ 1.2. Since that time scale is 
not explored by the Hinode data-set, we cannot exclude a 
terrestrial atmosphere effect: despite the MFBD and the 
jitter correction, the G-band bright points may undergo 
deformations and displacements as a result of still uncor- 


Table 1. 7 values retrieved by power law fits form the 
displacement spectra shown in Figs. |T] and [2] in different 
space and time ranges. 

Space Range [Mm] Time Range [sj 7 value 
< 0.05 <40 ~ 1.2 

0.05 - 1.50 40 - 1600 ~ 1.4 

> 1.50 > 1600s ~ 1.2 


reeled seeing aberrations. The seeing short (<C 1 s) corre¬ 
lation times could affect the results. On the other hand, 
the G-band bright points are not perfect passive tracers: 
they have a structure by their own, which can be deformed 
by the surrounding plasma flows. More, buffeting-induced 
transv erse oscillations as those found by IStanEfalini et al. 
( 20131) could generate a signal likely to be detected as a 
movement of the tracer. Those effect could mimick a ran¬ 
dom walk behaviour of the structure centre, for spatial 
scales < 100 km. The third point we raise is that the 
short time range t < 40 s is lower than the correlation 
time-scale of granulation, suggesting that, at least on that 
short scales, 7 ~ 1.2 can really be the signature of anoma¬ 
lous diffusion by the intergranular turbulent flow field. 
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What happens for scales longer than ~ 40 s and larger 
than ~ 100 kml Why 7 increases and indicates a clear 
super-diffusive regime? 

At those spatial and time scales all the points raised above 
do not hold any more, and the magnetic elements be¬ 
have like passive tracers on the large scale velocity field. 
Therefore, we interpret the 7 > 1 regime as a real signal 
of super-diffusivity. Usually, this is explicated as the re¬ 
sult of a turbulent dispersion process and analyzed in the 
framework of a fully developed turbulence. Instead, here 
we put forward the hypothesis that such an effect could 
be mainly due to the deformation of the lattice-like down¬ 
flow network. Therefore, the displacement spectrum would 
also embed the contribution of the deformation of the net¬ 
work, which could be dominant at the large scales. To test 
this hypothesis, we developed a simplified computational 
model to simulate such a regime. 

4. Simulation 

To test our hypothesis and investigate the possible con¬ 
tribution of an advecting granular field in diffusing the 
tracers, we developed a simple simulation and studied the 
motion of the tracers embedded in. 

Since the tracers used in the diffusion studies are often (if 
not always) embedded in the downflow network, we devel¬ 
oped an algorithm mimicking the advective flow generated 
by a granular field similar to the solar photosphere. In 
this simulation, the tracers are the vertexes of a Voronoi 
Tessellation generated from points representing the gran¬ 
ule upflow centres. Such vertexes are simply advected by 
the flow generated by the granules’ evolution, and no diffu¬ 
sive term is represented in the simulation. The interaction 
between the upffows of the granules is simulated by a re¬ 
pulsive force acting between the centres as a function of 
the associated upflow intensity. Stronger upffows tend to 
push away and ’squeeze’ weaker ones, causing an evolu¬ 
tion of the Voronoi tessellation, which in this simulation 
represent the lattice-like downflow network. To take into 
account the granules which are pushed outside of the sim¬ 
ulation box or are ’squezeed to death’ by their neighbours 
(a Sink term), we introduced a weak Source term as the 
generation in a random position of new upflows to main¬ 
tain the number of granules constant. The motion of the 
tracers therefore is defined only by the deformation of the 
lattice-like pattern specified by the upffows’ position and 
their relative intensity. 

The rules of the simulation are the following: 

1 ) At t = 0, V = 1024 points are randomly distribuited 
on a 2048 x 2048 matrix, such that the average distance 
between them is 64 pixels. 

2) Each point p has an associated interaction range R{p) 
expressed in pixels. At t = 0, we initialize randomly i?(p) 
with a Gaussian distribution centred at 64 pixels and with 
a standard deviation of 6 pixels. 

3) Every cycle, any point repels away the other points 
within its R{p) by 1 pixel. 

4) At the end of the cycle, each point has its R{p) de- 



d [arcsec] 

Figures. A snapshot from the simulation. The black 
crosses represent the advection centres; the red vectors 
are the displacement (5x exaggeration) to be applied to 
the centres due to the neighbours’ action; the red boxes on 
the Voronoi Tesselation highlight the cell vertexes, which 
are used as tracers to compute the displacement spectrum. 
For the sake of visualization, only about a twentieth of the 
simulation domain is shown. 


creased by the numbers of displacement it has undergone 
during the cycle due to the repulsion of the other points. 

5) If a point is displaced out of the simulation bounds or 
its R{p) falls below 1, it is replaced by another point in a 
random position and with R{p) = 64 ± 6 . 

6 ) At the beginning of any cycle, all the R{p) are incre¬ 
mented by a connectivity factor (7 = 4, to compensate for 
the displacements caused by 4 neighbours, on average. 

7) The algorithm is iterated for 1000 cycles. 

With the suitable choice of time-step {dt = 15 s) and 
pixel-scale (px = 1/64 arcsec)^ this algorithm simulates 
the interaction between the upflows of the solar granula¬ 
tion. The granules in a underdense location will expand, 
growing a larger and larger i?(p), while granules in a more 
crowded location may shrink and disappear, due to the 
other granules’ pressure. After some iterations, the simu¬ 
lation reproduces a statistically steady, renovating gran¬ 
ular pattern (a sample snapshot is shown in Fig. ED. We 
chose N, R{p)t=o = [64 ± 6 ] and C = 4 to minimize this 
transient phase and to obtain an average granular size of 
^ 1 arcsec = 750 km on the solar surface. The choice of 
such values, corresponds to assume that, on average, any 
granule interacts only with its first neighbours, since the 
average distance between centres is equal to the mean in¬ 
teraction range. 

After the transient phase, on average, 15 granule centres 
have to be generated randomly in every cycle to compen- 
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sate those granules which moved out of the computation 
bounds or were suppressed. Thus, the simulation contains 
a weak Sink/Source term, since the loss or the injection of 
a granule centre changes the network pattern locally and 
cancels out /generates new Voronoi vertexes and therefore 
new tracers. 

To be as consistent as possible with the real data analy¬ 
sis, we applied to the vertexes the same tracking algorithm 
used for the two data-set of Figs. [1] and [21 
We run several simulations, using different parameter 



10 100 1000 10000 
t [s] 

Figure 4. The displacement spectrum for the Voronoi 
Tessellation vertexes of the advection simulation. The lines 
represent power law fits for the ranges [t < 135 s] and 
[405 s<t< 1000 s]. 

sets, obtaining consistently super-diffusive 7 values, which 
are compatible with thos e retrieved by the observations 
( Giannattasio et al ][20ll). The displacement spectrum re¬ 
trieved by a sample simulation is shown in Fig.jT] The ex¬ 
ponent retrieved by a power law fit for times shorter than 
135 s is 7 ~ 1.5. For times longer than 400 s, the 7 value 
decrease to ~ 1.1. 

The transition point between the two behaviours has been 
estimated via a fitting procedure to be 200 ± 90s. The 
errors on the 7 values reported in the figure represent 
the variability of the fit results due to changes in the fit¬ 
ting ranges and not the fitting errors themselves, which 
are tipically an order of magnitude smaller. Changing the 


timestep dt or the pixel scale has no effect on the index 7. 
Changing instead the number of granules N or the con¬ 
nectivity C or the shape of the distribution of R{p) at 
t = 0, affects the average granule size and the duration 
of the transient phase. The 7 values themselves depend 
weakly on the parameters N, R{p) and C, which define 
the ratio between the average interaction range and the 
average granule dimension. The changes of the 7 values 
due to small variations of the parameter sets are included 
in the reported errors. Instead, those parameters affect 
the position of the transition between the 7 ~ 1.5 and 
7 ~ 1.1 regimes, since they set the spatio-temporal scale 
over which the flow can be considered isotropic: for a scale 
larger than average granule dimension and longer than the 
average granule duration. 

5. Simulation Discussion 

The motion of the vertexes of a lattice-like pattern, un¬ 
dergoing deformation due to the mutual repultion of the 
generating points, has a displacement spectrum with a 
spectral index 7 ~ 1.5. Given the parameters chosen for 
the simulation, the pattern produces a behaviour consis¬ 
tent with the solar granulation: it shows cells that expand, 
collapse and shuffle with characteristic time, velocity and 
dimension close to the granular ones. The displacement 
spectrum computed from the motion of the “downflow” 
vertexes of this simulation is totally defined by the pat¬ 
tern evolution, which has a single stochastic ingredient: 
the replacement of the removed upflows with new ones 
in random positions. Since we introduce only ~ 10 new 
upflows on the ^ 10^ existent upflows every step, we esti¬ 
mate this should be a small contribution. We recall that 
the number of new upflows introduced at each iteration is 
a result of the Voronoi dynamics of the simulation (not a 
chosen parameter), which is determined by the set of the 
parameters discussed above. Therefore, the 7 value ob¬ 
tained should derive only from the mutual interaction of 
the cells, which is nor entirely stochastic (this would lead 
to Brownian motion and 7 = 1), neither ballistic (7 = 2). 

6. Conclusions 


Several works in the literature utilized a Lagrangian ap¬ 
proach to determine the spectral index 7 in the displace¬ 
ment spectra (S'^) ^ of passive tracers on the solar 


photosphere (e.g. 

iRuzmaikin et al. 19961 Utz et al.l 20ld: 

Abramenko et al.l 

2011; Kitiashvili et al. 2012f). The most 

recent among; these works fsee Jafarzadeh et al. 20141 and 


references therein), agree in retrieving a 7 > 1 value, sug¬ 
gesting that a super-diffusion regime is present in the solar 
photosphere. 


In this paper, we computed with the same algorithm the 
displacement spectra of magnetic elements from two spec- 
tropolarimetric data-sets: one from a ground based instru¬ 
ment (IBIS at NSO/DST) and the other from a space- 
borne instrument (NFI on-board HINODE satellite). We 
confirm the superdiffusive motion of these elements, but 
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put forward an alternative interpretation in the term of 
the deformation of the lattice-like structure advecting the 
tracers, instead of a fully developed turbulence framework 
used so-far. 

To substantiate this interpretation, we realised a simple 
simulation of an advection process mimicking the granu¬ 
lar flows on the solar surface and used the vertexes of the 
simulated granular cell as tracers. In this simulation, the 
motion of the tracers is entirely determined by the defor¬ 
mation of the lattice-like downflow network modeling the 
solar surface flow patterns. When analysed with the same 
approach as the real data-sets, the simulation showed a 
displacement spectrum with a “super-diffusive” (7 ~ 1.5) 
behaviour. 

This result confirms that other theoretical frameworks 
may be used to interpret the motion of the magnetic ele¬ 
ments on the solar surface, and that the diffusion parame¬ 
ters which are calculated from the displacement spectrum 
analysis, may not be necessarily related to the effects of a 
turbulent reg ime. _ 

For example, Escande &: SattinI ( 2007ll investigated theo¬ 
retically some cases when the v term in the advection- 
diffusion equation cannot be zeroed by changing reference 
frame (the Lagrangian approach): in the hypothesis of an 
evolving advection field, the tracers can jump from one 
cell to the other and enter into long flights, depending on 
the characteristic temporal, spatial and intensity scales of 
the advection field. Such a point of view is definitely not 
far from t he simple simulation w e put forward. 

Parallely, ISolomon et al. (1994) explored experimentally 
different regimes, from quasi-periodic to chaotic and to 
turbulent regimes, finding that both the quasi-periodic 
and the chaotic regimes are characterised by Levy-flights 

and show super-diffusion, with 7 ~ 1.5.^ _ 

Other anal ogies between the vie ws of Escande fc SattinI 
( 2007ll and Solomon et al. (1994) and our algorithm and 
the possibility to further extend their r esults (possibly, ap¬ 


plying the pair separation approach of lOiannattasio et al. 


2014bfl . are at the moment under investigation. 
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